An Anderson-localized random nanolaser 
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Precision is a virtue throughout science in general and in optics in particular where carefully fab- 
ricated nanometer-scale devices hold great promise for both classical and quantum photonics In 
such nanostructures, unavoidable imperfections often impose severe performance limits but, in cer- 
tain cases, disorder may enable new functionalities/ Here we demonstrate on-chip random nanolasers 
where the cavity feedback is provided by the intrinsic disorder in a semiconductor photonic-crystal 
waveguide, leading to Anderson localization of light P This enables highly efficient and broadband 
tunable lasers with very small mode volumes. We observe an intriguing interplay between gain, 
dispersion-controlled slow light, and disorder, which determines the cross-over from ballistic trans- 
port to Anderson localization. Such a behavior is a unique feature of non-conservative random media 
that enables the demonstration of all-optical control of random lasing. Our statistical analysis shows 
a way towards ultimate thresholdless random nanolasers. 



The last decade has witnessed important breakthroughs 
in the search for the ultimate nanolaser with zero threshold 
and ultra-small size^ ^ . This quest involves perfecting the 
nanocavity that mediates the optical feedback, but any fabri- 
cation method has finite accuracy making disorder ubiquitous, 
which induces light loss and results in hampered functionali- 
ties. Random lasers may overcome this limitation since they 
originate from the random multiple scattering of light induced 
by disordei^^^'^. However, their multidirectional lasing emis- 
sion, the lack of a broadband emission tunability, and the 
poor lasing mode confinement found in standard - diffusive 
- random lasers give rise to strong mode competitiorP^ and 
chaotic laser emissiorP', and has so far limited their prac- 
tical usability. Here we report on broadband-tunable ran- 
dom lasing in photonic-crystal waveguides in the regime of 
Anderson localization with low thresholds and small mode 
volumes. The strong confinement achieved by Anderson lo- 
calization reduces the spatial overlap between lasing modes, 
thus preventing mode competition, which improves the laser 
output stability. Furthermore, a rich physical phase space 
is unraveled where the system can be driven in and out of 
the Anderson- localization regime by introducing gain or loss, 
thus demonstrating all-optical control of the lasing process. 
A statistical analysis of the data reveals the pathway towards 
minimizing the lasing threshold and mode volume by reducing 
the localization length, which provides unique possibilities of 
controlling and optimizing random lasing at the nanoscale. 

The intrinsic disorder present in state-of-the-art photonic- 
crystal waveguides turns out to be sufficient to induce 
Anderson-localized cavities in the slow-light regime of waveg- 
uide propagation. Figure [TJl illustrates such spontaneous las- 
ing cavities that are spatially distributed along the waveg- 
uide and spectrally bound to the slow-light regime^^. In 
the present experiment, quantum wells are embedded as gain 
medium in the photonic crystal membrane and are optically 
pumped. For low pump power, the strong absorption in the 
quantum wells damps interference and hence also Anderson 
localization, as illustrated in Fig. [TJ), while for increasing ex- 
citation power, light amplification compensates the losses and 
stimulates the transition to the random lasing regime, as il- 
lustrated in Fig.[l]:. Figure[T]l displays the emission spectrum 
of a single Anderson-localized lasing mode below and above 



threshold, while Fig.[T^ shows the peak output emission inten- 
sity vs. input excitation power. A distinct lasing threshold is 
observed, accompanied by a decrease of the cavity line width 
with excitation power that are the signatures of laser oscilla- 
tion. Two important figures-of-merit of a nanolaser are the 
/3-factor and the cavity mode volume, which are the fraction 
of spontaneous emission that couples to the lasing mode and 
the spatial extension of the lasing mode, respectively. They 
can be determined by modeling the laser input-output curves 
with micro-cavity semiconductor-laser rate equations^^ (see 
Fig. [l^). Details on the model can be found in the supple- 
mentary information (S.I.). For this particular lasing mode 
we obtain a /3-factor of 0.31, which interestingly is signifi- 
cantly larger than previous values reported in photonic-crystal 
lasers with quantum- well gain media^*^. A mode volume of 
4.53 (A/n)^ is obtained, where A = 1587 nm is the lasing wave- 
length and n = 3.4 is the refractive index of InGaAsP, which is 
consistent with the mode extent found for Anderson- lo calize d 
modes in cavity quantum electrodynamics experimentP^^. 

Determined by the size of the pumped area on the waveg- 
uide, multi-frequency lasing can be observed, as presented in 
Fig. [if and[l^. As shown in Fig. [if , the onset of multi-mode 
lasing tends to be sequential, i.e., with increasing excitation 
power the laser peaks grow one by one, each displaying the 
characteristic input-output lasing curve (see Fig.[l^). No sig- 
nature of mode competition is observed, which is attributed to 
the fact that the Anderson-localized lasing modes tend to be 
so strongly confined that different modes do not deplete the 
same gain region. In diffusive random lasers, on the contrary, 
the extended lasing modes imply that different modes overlap 
spatially, which induces strong mode competition and chaotic 
lasing emission^^. The improved stability of the multi-mode 
Anderson-localized random laser compared to diffusive ran- 
dom lasers may have significant advantages for mode-locked 
operation^^ . 

In diffusive random lasers, the lasing wavelength is primar- 
ily determined by the spectral peak of the gain mediunP^ 
while only weak dispersion has been incorporated into the 
multiple-scattering medium enabling modest tunability of the 
random lasei^^ A photonic-crystal waveguide, on the con- 
trary, is highly dispersive, which implies that the Anderson- 
localized random lasers can be spectrally tuned over an un- 
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Figure 1: Single-mode and multi-mode random las- 
ing in the Anderson-localization regime, a, Illustration 
of a random laser in a disordered photonic-crystal waveguide 
formed from Anderson- localized modes, where the different 
colors indicate that multifrequency lasing can occur, b, and 
c, Illustration of multiple-scattering process below and above 
threshold, respectively. Below threshold, loss is dominat- 
ing and Anderson-localized cavities do not build up. Above 
threshold, gain compensates loss and Anderson- localized ran- 
dom lasing sets in. d. Measured photoluminescence spec- 
tra of a single Anderson-localized lasing mode below (lower 
plot) and above (upper plot) the lasing threshold, e. Output 
intensity (black squares) and linewidth (blue squares) mea- 
sured vs. excitation power. The solid line represents a fit to 
semiconductor-laser rate equations and the black dots repre- 
sent the background quantum- well fluorescence displaying no 
lasing. The arrows indicate the excitation power of the spec- 
tra shown in d. f, Photoluminescence spectra vs. excitation 
power in multi-mode operation where the different colored 
curves are fits with Lorentzian functions representing differ- 
ent lasing modes, g. Output intensity vs. excitation power 
corresponding to the different lasing modes shown in f. The 
solid lines represent the fit to semiconductor-laser rate equa- 
tions for each lasing mode. 



precedented wavelength range by varying the photonic-crystal 
lattice, with a tunability that is only limited by the band- 
width of the quantum well gain medium (almost 300 nm). 
Figure [2^ shows the spectral tuning of the lasing modes over 
50 nm when varying the ratio r/a, where r is the hole radius 
and a the lattice constant of the photonic-crystal waveguide. 
This broad tunability is due to a particular characteristic 
of Anderson-localized modes in photonic-crystal waveguides: 
the modes appear solely in the slow-light regime of the waveg- 
uide near the band-edge cutoff, where the density of optical 
states (DOS) associated with the waveguide mode is large^^. 
Figure [2]3 displays how the high-DOS region of the waveguide 
can be tuned by varying the photonic-crystal parameters, es- 
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Figure 2: Tunable Anderson-localized random lasing. 

a, Photoluminescence spectra (vertically shifted for visual 
clarity) above lasing threshold collected for various values of 
r/ a, where r is the hole radius and a the lattice constant of 
the photonic crystal. The quantum well emission spectrum is 
shown for comparison (black line and inset), b. Calculated 
density of states of a photonic-crystal waveguide without dis- 
order where the different colors correspond to the different 
values of r/ a shown in a. 



tablishing the tuning mechanism for the Anderson-localized 
modes. 

Anderson localization is a random process and is described 
only by a statistical distribution of lasing parameters. The 
distribution of Q-factors at lasing threshold for all the modes 
observed in the experiment is shown in Fig. [3^, from which 
we can extract the localization length and the loss lengt 
The localization length is the average distance between each 
random scattering event and we extract ^ = 6.0 ^im together 
with a loss length of / = 800 \im. The latter is determined 
by out-of-plane scattering and absorption in the unpumped 
region of the quantum well, and is likely to be biased towards 
low loss since preferentially the very localized modes will lase. 
Since the localization length is significantly shorter than the 
sample length of 100 \im, the sample is found to be deeply 
in the Anderson- localization regime at lasing threshold. Fig- 
ure [Sj) displays the distribution of lasing thresholds observed 
experimentally. An asymmetric threshold distribution with 
the long tail extending to high threshold values is obtained as 
predicted by theory^, in contrast to the symmetric threshold 
distributions observed for diffusive random lasers^^. The av- 
erage threshold observed of 418 \iW is comparable to the typ- 
ical threshold of standard photonic-crystal lasers pumped in 
a similar geometr}^^^, illustrating that Anderson-localized 
cavities may compete favorably with engineered nanocavity 
lasers that are hampered by imperfections. The /3-factor and 
mode- volume distributions are plotted in Fig. [3]: and Fig.jsjl. 
Large /3-factors ranging from 0.28 to 0.38 are found for all las- 
ing modes in conjunction with small mode volumes ranging 
from 4 to 6 (X/n)^ . The mode- volume distribution enables an- 
other estimate of the mode-length extension along the waveg- 
uide (see S.I.), which is plotted as the top horizontal axis in 
Fig. [Sji. The average mode length of 4.3 ^im is very similar 
to the localization length (^ = 6.0 \im) that is extracted from 
modelling the Q-factor distributions. This remarkable agree- 
ment confirms the validity of the two independent models 
applied for analyzing the experimental data, and illustrates 



3 



Mode length (|im) 

4.0 4.5 



5.0 






2000 4000 6000 8000 

Q - factor 



200 300 400 500 600 700 

Threshold (^W) 





700 -| 




600- 






■D 


500- 


O 






400- 






1- 


300- 




200 




0.27 0.30 0.33 0.36 

(5- factor 

0.39 



0.39 



4.5 5.0 5.5 

Mode volume (Vnf 




2000 4000 6000 8000 

O - factor 



Mode volume (^nf 



2000 4000 6000 8000 

O - factor 



4.5 5.0 5.5 

Mode volume (X/n) 



Figure 3: Statistical properties of random lasing in the Anderson-localization regime, a, Experimental Q-factor 
distribution of the observed localized modes in the photonic-crystal waveguide at the threshold power (histograms) fitted with 
a ID multiple-scattering model to extract the localization length (black line), b, Laser-threshold distribution obtained from 
measured input-output curves, c, /3-factor and d, mode-volume (converted to mode length on the top axis) distributions 
extracted from measured input-output curves. e,f and g,h Measured threshold and /3-factor versus Q-factor and mode volume, 
respectively. The Spearman correlation coefficient, rg, describes the correlation of the different parameters, with 1 corresponding 
to perfect correlation and -1 to perfect anti-correlation. The dashed lines are guides to the eye. 



that important microscopic parameters of Anderson-localized 
random lasers can be reliably extracted. 

An important question is how to improve the lasing figures- 
of-merit in the case of random lasing. The ultimate nanolaser 
will ideally be thresholdless due to the /3-factor approaching 
unity. The statistical analysis presented here offers insight 
to the governing parameters. Figures [3^,f and |3^,h show 
the measured lasing threshold and /3-factor versus Q-factor 
and mode volume. The correlation between these parameters 
is quantified by the Spearman's rank correlation coefficient 
rg (details in S.I.) that demonstrates how a high Q-factor 
and a small mode volume improve lasing performance. In 
Anderson-localized random lasers these cavity characteristics 
are directly linked to the localization length scaled to the 
sample length, which is a universal parameter determining 
all properties of light transport in the random medium. The 
key to improving the lasing performance of the Anderson- 
localized lasers even further would be to shorten the localiza- 
tion length, which may be feasible, e.g., by combining ran- 
dom disorder with long-range correlated order. Indeed, a 
shorter localization length would be beneficial in two ways 
since it would simultaneously increase the laser Q-factor and 
decrease the mode volume^. However, unlike transport prop- 
erties of random media, the Anderson-localized random lasers 
are not fully determined by universal parameters, and the 
non- universality of light emission in a random medium has re- 
cently been shown experimentall}^. The current exploration 
of light emission and propagation in random media reveals 



the exciting prospects of disordered nanostructures for lasers 
and cavity quantum electrodynamics experiments. 

Any amount of absorption or light leakage from the struc- 
ture reduces the confinement due to Anderson localization. 
While inclusion of large amounts of gain inevitably modifies 
the properties of the underlying ordered structure^', moder- 
ate optical gain can compensate the effect of loss, thereby 
recovering the confinement of light. The effect of gain on An- 
derson localization can be investigated by extracting the prob- 
ability distribution of the intensity fiuctuations P(/ = //(/)), 
where (/) denotes the intensity obtained after ensemble av- 
eraging over configurations of disorder. The intensity, /, 
is measured while raster scanning the sample, as shown in 
Fig. [4^ (see S.I. for details), and displayed for different exci- 
tation powers in Fig.|4]3. For low excitation power, the strong 
absorption from the quantum wells damps long interference 
paths resulting in the normal distribution of P(/) shown in 
Fig. |4l3. When increasing the excitation power, light am- 
plification compensates absorption, leading to Anderson lo- 
calization and random lasing, for which P(/) displays the 
characteristic log-normal distribution of Anderson localiza- 
tion. Although no complete theory of Anderson localization 
in active media is yet available, the cross-over to localization 
can be qualitatively reproduced with a ID multiple scattering 
model^^ where the effect of gain is modeled by increasing the 
loss length, see Fig. |4]:. The experimentally recorded vari- 
ance of the intensity fiuctuations, Var(/), is found to increase 
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crystal waveguide, tailored to the non- localizing "fast light" 
regime, next to the Anderson-localizing waveguide, which 
would enable a highly efficient and stable multi-color laser 
with directional output. 
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Figure 4: All-optical control of the cross-over to An- 
derson localization, a, Photoluminescence intensity for 
an excitation power of 500 ^lW collected while scanning the 
excitation and collection objective along a photonic-crystal 
waveguide, b. Intensity probability distribution for different 
excitation powers, c. Calculated intensity probability distri- 
bution with the model of Ref. |21| with a localization length 
of ^ = 6 \im and two different loss lengths. / = 800 ^im corre- 
sponds to the loss length extracted from the Q-factor distri- 
butions measured at threshold, d. Variance of the normalized 
intensity as a function of the excitation power. The horizontal 
dashed line shows the localization criterion of Var(/) = 0.5, 
and the vertical line indicates the excitation power at which 
the cavities start lasing. 



monotonically with excitation power, as shown in Fig.[4]l. At 
the onset of random lasing, Var(/) = 0.85 is recorded, which 
is larger than the threshold for Anderson localization in a 
ID structurd^^ of Var(/) = 0.5. The ability to manipu- 
late loss and gain in a non-conservative ID random system 
opens for the possibility of exploring the very rich parame- 
ter phase space that has been theoretically proposed^"^. The 
all-optical control of the cross-over from a lossy regime to An- 
derson localization offers exciting possibilities of fine-tuning 
the Anderson-localized modes. 

Embedding light emitters in Anderson-localized random 
media provide a promising new pathway to enhanced light- 
matter interaction relying on the natural occurrence of cav- 
ities in random media rather than nano-engineering. The 
rich and intricate physics is currently being explored, but 
the prospective applications for low-threshold nanolasers and 
highly- efficient single-photon sources appear very attractive 
since the relevant figures-of- merit for these applications match 
engineered systems, with the added value that random struc- 
tures could be much easier to fabricate. The fundamental lim- 
its of the disordered media are still to be explored, but the per- 
formance is likely to be improved even further by deliberately 
mixing the amount of order and disorder in the structures. 
An important issue for applications of an Anderson-localized 
nanolaser will be the ability to couple light from the randomly 
positioned laser to a well-defined single mode. This could eas- 
ily be implemented on-chip by positioning another photonic- 
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SUPPLEMENTARY INFORMATION 

Sample fabrication 

The semiconductor heterostructure used for the fabrication 
of the photonic-crystal waveguides was grown on a (100)- 
oriented semi-insulating indium phosphide (InP) substrate 
by metalorganic chemical vapor deposition (MOCVD). The 
structure consists of a 340-nm-thick Ino.77GaAso.503P mem- 
brane, which is grown on a 100-nm-thick InP, 100-nm-thick 
InAlAs and 800-nm-thick InP sacrificial layer, incorporating 
10 layers of Ino.75GaAso.86P quantum well at the center with a 
thickness of 5 nm and a separation of 5 nm between adjacent 
quantum wells. The total refractive index of the structure 
is n = 3.4. The photonic-crystal waveguides are fabricated 
by electron-beam lithography, reactive-ion etching, and wet 
etching using hydrochloric and hydrofluoride acid solutions. 
Samples with lattice constant a — 380 nm, and a range of dif- 
ferent hole radii 0.237a < r < 0.263a are fabricated. Finally, 
340-nm-thick suspended membranes are formed by removing 
the sacrificial layer. 



Experimental setup and optical characterization 

Microphotoluminescence measurements are carried out at 
room temperature. A 120 fs pulsed Ti:sapphire laser oper- 
ated at 800 nm with a repetition rate of 80 MHz is used for 
optical pumping. An excitation beam is focused to a spot 
with a diameter of 1.5 |J,m on the surface of the sample using 
a 50x microscope objective lens (numerical aperture: 0.65), 
and aligned to the photonic-crystal waveguides using piezo- 
electric nanopositioners. The photoluminescence signal is col- 
lected by the same objective lens within a diffraction-limited 
region in a wide wavelength range 1500 nm < A < 1600 nm, 
dispersed by a 300 mm grating spectrograph with a spectral 
resolution of 0.1 nm, and detected using a liquid-nitrogen- 
cooled InGaAs CCD (charge-coupled device) camera. 

We measure a total of 25 input-output curves by vary- 
ing the excitation and collection position along the photonic- 
crystal waveguide. The laser threshold is obtained by ex- 
trapolating each input-output curve in a linear scale to zero 
power. The Q-factor is extracted by fitting the cavity mode 
with a Lorentzian at the threshold power. Two out of the 
25 input-output curves show the cavity mode also below the 
laser threshold (Fig. [T]3 shows one example). For the rest of 
the measurements, the excitation power at which the cavity 
mode appears matches the threshold power extracted from 
the input-output curve fitting. Based on this fact, we can 
increase the statistics of threshold and Q-factor distributions 
by gradually increasing the excitation power and recording 
the power at which the mode appears and its corresponding 
spectrum. By doing so, 25 extra data points for the laser 
threshold and Q-factor are added without recording the full 
input-output curve. The full data set is plotted in Fig. [3]: and 
111. 

The intensity probability distribution P{I) is measured 
by collecting the intensity I^^y while scanning the excitation 
and collection objective along the photonic-crystal waveguide 
at each spatial position (x,^) with a spatial binning size of 
0.25 ^Lm and at different wavelengths with a binning size of 
1 nm. Subsequently, an average over the wavelength range 
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1580 nm — 1595 nm is performed. 



Rate-equation analysis 

We use a modified semiconductor-laser rate-equation model 
for the carrier density N and the photon density P of the 
random lasei^^ 



dN 



— Rry 



dt 

dP on P 

^ = ri3BN^ + r^G{N)P , 

dt rp 



G{N)P -AN - CN^ 



(1) 
(2) 



where Rp is the pulsed injected carrier density with the tem- 
poral width of the excitation laser pulse, B is the bimolec- 
ular recombination coefficient of the unstructured material 
that we assume here to be unchanged by the photonic crys- 
tal structure, and A and C are the non-radiative coefficients 
corresponding to surface recombination and Auger processes, 
respectively F and are the conventional and generalized 
confinement factors, G{N) — c/nef£Go{N-Ntr) is the gain co- 
efficient, and Tp is the photon lifetime obtained from the mea- 
sured Q-factor at the laser threshold (the remaining parame- 
ters are defined in Table [l]). The formalism takes into account 
Purcell enhancement of the spontaneous and stimulated emis- 
sion^^ including the fact that the quantum well electronic 
density of states is very broadband compared to the cavity 
bandwidth^^. Furthermore, in the present experiment where 
only a small volume is optically pumped, the cavity mode vol- 
ume is comparable to the volume of the gain material. The 
conventional optical confinement factor, F = Va/V, where 14 
and V are the active material and modes volume, respectively, 
does not take into account the spatial-modal overlap profile 
with the gain, and a generalized confinement factoi^^ should 
be used. The mode volume V and the generalized confinement 
factor F*^ are defined as 



V 



/e.(r)|Ee(r)|^d^r 
er(ro)|Ee(ro)|2 
/^^6.(r)|Ee(r)|^d^r 
- /er(r)|Ee(r)Pd3r 



(3) 
(4) 



where ro is the position of the anti-node of the cavity field, V 
is the mode volume of the Anderson-localized laser mode, Va 
is the active volume covered by the carriers after excitation 
and subsequent diffusion, and Ec(r) is the cavity-mode elec- 
tric profile. In order to reduce the number of free parameters, 
we simplify the expression of F^ by assumine: that the cav- 
ity field is uniformly distributed spatially with the amplitude 
|Ec(r)|2=||Ec(ro)|^ and er(r) = er(ro), leading to 



G 



V 

^ a 

2V 



(5) 



where VCp is the overlap volume between the cavity mode and 
the active volume. 

We estimate this overlap volume as follows. The size of 
the active area is determined by the excitation beam, with a 
1.5 ^im diameter (red area in Fig. [5^). The carrier diffusion 
distance is / = V Dr = 0.9 ^im, by assuming D — 8.0cm^/s 
as the carrier diffusion coefficient and r — 1 ns as the car- 
rier lifetime. We neglect the carrier diffusion distance in the 



mVs 
X 10-^^ mVs 



TABLE L Material parameters for the semiconductor-laser 
rate-equation^^ 

Surface recombination rate {A) = 2.5 x 10^ s""*^ 
Bimolecular recombination rate [B) — 10 
Auger non-radiative recombination rate (C) = ( 

Transparency carrier density (A^tr) = 10^^ m 
Gain coefficient (Go) = 2.5 x 10^ m^ 
Effective refractive index (neff)= 2.8 



estimation of the active area since the carrier lifetime is ex- 
pected to be much shorter above threshold giving rise to a 
much shorter diffusion distance. We assume the Anderson- 
localized modes to be confined transversally by the effect of 
the 2D photonic-crystal gap (blue area in Fig.|5^). Therefore, 
a total overlap length of 1.5 ^im is determined by the active 
area, while the overlap width of 2a sin ^ — 2r — 458 nm is es- 
timated by the width of the photonic-crystal waveguide. The 
overlap height is given by the total thickness of 10 QW layers 
(5nm X 10). In summary, we estimate 



(1500 X 458 X 50) 



14 = ((tt X 750^ - 14 X TT X 100^) x 50) nm^ 



(6) 
(7) 




Figure 5: Schematic illustration of the active area and 
cavity mode area, a. Scanning electron micrograph of a 
photonic-crystal waveguide. The red area indicates the op- 
tically excited and thus active area, 14, and the blue area 
represents the area of a random laser mode, determining V . 
The overlap between them gives the effective active over- 
lap, . b, 2D finite- difference time-domain calculations of 
the ^-component of the electromagnetic field in a disordered 
photonic-crystal waveguide with rj a — 0.28 and b — 3%. 
Strongly confined Anderson-localized modes are found in the 
slow-light regime after exciting the structure with a broad- 
band source centered at 0.26a/A with a bandwidth of 0.02a/A. 
The intensity of the ^-polarized electric field is plotted on top 
of the simulated structure. 
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where the area of 14 holes is subtracted from the excitation 
spot area to account for the final active area (see Fig. |5^). 

Finally, we estimate the cavity mode length by dividing 
the mode volume obtained from the semiconductor-laser rate- 
equation model with an effective transverse mode area Aeff. 
To calculate Aeff , we assume the Anderson-localized modes to 
be confined transversally. The width of the cavity mode can 
be calculated from finite-difference time-domain simulations 
of a photonic-crystal waveguide where disorder is induced by 
randomly varying the hole positions according to a Gaussian 
distribution with a standard deviation S = 3% relative to 
the lattice spacing (see Fig. ^p) and the thickness of the cav- 
ity mode is assumed to be the same as the membrane. By 
averaging over one unit cell along the waveguide we obtain 
Aeff - 0.11 [im^. 



Investigating correlations of random laser 
parameters 

In order to investigate how to optimize the performance 
of a random laser, we estimate the Spearman's rank correla- 
tion coefficient, rs, between the different laser parameters. In 



particular, we focus on the dependence of the laser thresh- 
old Pth and the /3-factor on the Anderson-localized mode 
properties, i.e., the Q- factor and mode volume, rs is a non- 
parametric measure of the statistical correlation between two 
variables, regardless the particular functional dependence be- 
tween them^^. The Spearman correlation varies from 1 (fully 
correlated) to -1 (fully anti-correlated) . 



* Present address: Universitat Konstanz, Fachbereich 

Physik, Fach M 700, 78457 Konstanz, Germany. 
^ Electronic address: ,lodahl@nbi.ku.dk^ URL: |www . | 
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